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■ The QCD phase transition at finite temperature is studied with the dual Ginzburg-Landau 

theory, which is the QCD effective theory based on the dual Higgs mechanism by QCD-monopole 
condensation. At high temperature, the confinement force is largely reduced by thermal effects, 
which leads to the swelling of hadrons. Simple formulae for the surface tension and the thickness 
| of the phase boundary are derived from the shape of the effective potential at the critical tem- 

perature. We investigate also the process of the hadron-bubble formation in the early Universe. 
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\Q 1. Introduction 

^Cj ■ Quantum chromodynamics (QCD) is the fundamental theory of the strong 

interaction [1-3]. In spite of the simple form of the QCD lagrangian, 
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<D . | 
^ ' £ Q CD = -rtlG^GT + q(P- m q )q, (1.1) 
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^ ■ it miraculously provides quite various phenomena like color confinement, dynamical 
chiral-symmetry breaking, non-trivial topologies, quantum anomalies and so on, as 
shown in Fig.l. Then, it would be interesting to compare QCD with the history of 
the Universe, because a quite simple 'big bang' also created various things including 
galaxies, stars, lives and thinking reeds [4]. Therefore, QCD can be regarded as 
an interesting miniature of the history of the Universe. This is the most attractive 
point of the QCD physics. 

Since it is quite difficult to understand the various QCD phenomena and their 
underlying mechanism at the same time, many methods and models have been 
proposed to understand each phenomenon. We show in Fig. 2 a brief sketch on the 
history of QCD and typical QCD effective models [1] as listed in Table 1. The 
bag model and the Nambu-Jona-Lasinio (NJL) model are simple relativistic quark 
theories. The nonlinear a model and the Skyrme-Witten soliton model are de- 
scribed by the Nambu-Goldstone pion. In particular, the Skyrme-Witten model 
is stimulative because a fermionic baryon can be described as a soliton made by 
bosonic pions in this framework. The origin of this magic is found in the non-trivial 
topology of the flavor dynamics in QCD. The lattice QCD simulation [5] using a 
supercomputer is a hopeful and promising method based on QCD directly, and 
its importance becomes larger and larger according to the great progress of the 
computational power. In particular, recent lattice QCD studies [6] shed a light 



on the confinement mechanism, which is one of the most difficult problems in the 
particle physics. In these years, the origin of color confinement can be recognized 
as the dual Higgs mechanism by monopole condensation, and the nonperturbative 
QCD vacuum is regarded as the dual superconductor. The dual Ginzburg-Landau 
(DGL) theory [7-14] was formulated with this picture, and can reproduce confine- 
ment properties like the string tension and the hadron flux-tube formation. 

Now, you may find a current of the QCD physics. In '80s, chiral symmetry 
breaking was the central issue. The chiral bag model, the NJL model and the a 
model were formulated with referring chiral symmetry. In '90s, on the other hand, 
the confinement physics is providing an important current of the hadron physics. 
The key word for the understanding of confinement is the "duality", which is re- 
cently paid attention by many theoretical particle physicists after Seiberg-Witten's 
discovery on the essential role of monopole condensation for the confinement in a 
supersymmetric version of QCD [15]. 

2. Color Confinement and Dual Higgs Mechanism 

We briefly show the modern current of the confinement physics. About 20 
years ago, Nambu, 't Hooft and Mandelstam proposed an interesting picture for 
color confinement based on the analogy between the superconductor and the QCD 
vacuum [16]. In the superconductor, the magnetic field is excluded due to the 
Meissner effect, which is caused by Cooper-pair condensation. As the result, the 
magnetic flux is squeezed like the Abrikosov vortex. On the other hand, the color- 
electric flux is excluded in the QCD vacuum, and therefore the squeezed color-flux 
tube is formed between color sources. Thus, these two systems are quite similar 
and can be regarded as the dual version each other. This idea is based on the 
"duality" of the gauge theories, which was firstly pointed out by Dirac more than 
50 years ago [17]. 

With referring Table 2 and Fig.3, we compare the ordinary electromagnetic 
system, the superconductor and the nonperturbative QCD vacuum regarded as 
the dual superconductor. In the ordinary electromagnetism in the Coulomb phase, 
both electric flux and magnetic flux are conserved, respectively. The electric- 
flux conservation is guarantied by the ordinary gauge symmetry. On the other 
hand, the magnetic-flux conservation is originated from the dual gauge symmetry 
[8-14], which is the generalized version of the Bianchi identity. As for the inter- 
charge potential in the Coulomb phase, both electric and magnetic potentials are 
Coulomb-type. 

The superconductor in the Higgs phase is characterized by electric-charge con- 
densation, which leads to the Higgs mechanism or spontaneous breaking of the 
ordinary gauge symmetry, and therefore the electric flux is no more conserved. 
In such a system obeying the London equation, the electric inter-charge potential 
becomes short-range Yukawa-type similarly in the electro-weak unified theory. On 
the other hand, the dual gauge symmetry is not broken, so that the magnetic flux 
is conserved, but is squeezed like a one-dimensional flux tube due to the Meissner 



effect. As the result, the magnetic inter-charge potential becomes linearly rising 
like a condenser. 

The nonperturbative QCD vacuum regarded as the dual Higgs phase is charac- 
terized by color-magnetic monopole condensation, and resembles the dual version 
of the superconductor, where the "dual version" means the interchange between 
the electric and magnetic sectors. Monopole condensation leads to the spontaneous 
breaking of the dual gauge symmetry, so that color-magnetic flux is not conserved, 
and the magnetic inter-change potential becomes short-range Yukawa-type. Note 
that the ordinary gauge symmetry is not broken by such monopole condensation. 
Therefore, color-electric flux is conserved, but is squeezed like a one-dimensional 
flux-tube or a string as a result of the dual Meissner effect. Thus, the hadron flux- 
tube is formed in the monopole-condensed QCD vacuum, and the electric inter- 
charge potential becomes linearly rising, which confines the color-electric charges 
as quarks [8-11]. 

As a remarkable fact in the duality physics, these are two "see-saw relations" 
between in the electric and magnetic sectors. 

(1) There appears the Dirac condition eg = 2ir [17] in QED or eg = 4n [8] in QCD. 
Here, unit electric charge e is the gauge coupling constant, and unit magnetic 
charge g is the dual gauge coupling constant. Therefore, a strong-coupling system 
in one sector corresponds to a weak-coupling system in the other sector. 

(2) As shown in Fig.3, the long-range confinement system in one sector corresponds 
to a short-range (Yukawa-type) interaction system in the other sector. 

Let us consider usefulness of the latter "see-saw relation". One faces highly 
non-local properties among the color-electric charges in the QCD vacuum because 
of the long-range linear confinement potential. Then, the direct formulation among 
the electric-charged variables would be difficult due to the non-locality, which 
seems to be a destiny in the long-distance confinement physics. However, one 
finds a short-range Yukawa potential in the magnetic sector, so that the electric- 
confinement system can be approximated by a local formulation among magnetic- 
charged variables. Thus, the confinement system, which seems highly non-local, 
can be described by a short-range interaction theory using the dual variables, which 
is the most attractive point in the dual Higgs theory. 

Color-magnetic monopole condensation is necessary for color confinement in 
the dual Higgs theory. As for the appearance of color-magnetic monopoles in QCD, 
't Hooft proposed an interesting idea of the abelian gauge fixing [18,19], which is 
defined by the diagonalization of a gauge-dependent variable. In this gauge, QCD 
is reduced into an abelian gauge theory with the color-magnetic monopole [8,18], 
which will be called as QCD-monopoles hereafter. Similar to the 't Hooft-Polyakov 
monopole [2] in the Grand Unified theory (GUT), the QCD- monopole appears 
from a hedgehog configuration corresponding to the non-trivial homotopy group 
7r 2 (SU(A^ c )/U(l) iVc - 1 ) = Z^- 1 on the nonabelian manifold. 

Many recent studies based on the lattice gauge theory have supported QCD- 
monopole condensation and abelian dominance [6,13,20-22], which means that only 
abelian variables are relevant for the nonperturbative QCD, in the 't Hooft abelian 



gauge. Hence, the dual-superconductor scenario seems workable for color confine- 
ment in the QCD vacuum, and the nonperturbative QCD would be described by 
the dual Ginzburg-Landau (DGL) theory, which is the QCD effective theory based 
on the dual Higgs mechanism. 



3. Dual Ginzburg-Landau Theory 

The dual Ginzburg-Landau (DGL) lagrangian [13,14] for the pure-gauge system 
is described with the dual gauge field and the QCD-monopole field Xi 

Cdgl = tr {-^OV^ - dvBrf + [Dp, X ]t[^, X ] ~ Kx ] X - ^ 2 ) 2 } , (3.1) 

where = <9 M + igB^ is the dual covariant derivative including the dual gauge 
coupling constant g = 4n/e. 

The dual gauge field B^ = B^ ■ H = B 3 T 3 + B®T 8 is defined on the dual 
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gauge manifold U(l) m x U(l) m [13,14], which is the dual space of the maximal 
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torus subgroup U(l) e x U(l) e embedded in the original gauge group SU(3) C . The 
abelian field strength tensor is written as — *(d A B)^ so that the role of the 
electric and magnetic fields are interchanged in comparison with the ordinary 
description. 

The QCD-monopole field x i s defined as x = Yla=i XaE a [13,14] where 
E x = ^(T 6 + iT 7 ), E 2 = ^(T 4 - iT 5 ) and E 3 = ^(Ti + iT 2 ). Here, X a has 
the magnetic charge ga proportional to the root vector a. In the QCD-monopole 
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condensed vacuum with \xa\ = v, the dual gauge symmetry U(l) m x U(l) m is 

3 8 

spontaneously broken instead of the gauge symmetry U(l) e x U(l) e . Through 
the dual Higgs mechanism, the dual gauge field B^ acquires its mass = V%gv, 
whose inverse provides the radius of the hadron flux tube [8] , and the dual Meissner 
effect causes the color-electric field excluded from the QCD vacuum, which leads 
to color confinement. The QCD-monopole fluctuations Xa = Xa — v (a=l,2,3) also 
acquire their mass m x = 2y/\v in the QCD-monopole condensed vacuum. As a 
relevant prediction, only one QCD-monopole fluctuation x= Ea=l Xa appears as 
a color-singlet scalar glueball in the confinement phase, although the dual gauge 
field B^ and the other two combinations of the QCD-monopole fluctuation are not 
color-singlet and cannot be observed [13,21,22]. 

The DGL theory reproduces confinement properties like the inter-quark po- 
tential and the hadron flux-tube formation. We studied effects of the flux-tube 
breaking by the light quark-pair creation in the DGL theory, and derived the in- 
frared screened inter-quark potential [8-10], which is observed in the lattice QCD 
with dynamical quarks [5] . We studied also the dynamical chiral-symmetry break- 
ing (D^SB), which is also an important feature in the nonperturbative QCD, by 
solving the Schwinger-Dyson equation for the dynamical quark [8-10]. The quark- 
mass generation is brought by QCD-monopole condensation, which suggests the 



close relation between D^SB and color confinement. Thus, the DGL theory pro- 
vides not only the confinement properties but also D%SB and its related quantities 
like the constituent quark mass, the chiral condensate and the pion decay constant 
[8-10,23,24]. 

4. QCD Phase Transition in the DGL Theory 

In this chapter, we study the QCD phase transition [10-12,14,24] in the DGL 
theory. Although quarks and gluons are confined inside hadrons in the nonpertur- 
bative QCD vacuum, these colored particles are liberated and the system becomes 
the quark-gluon-plasma (QGP) phase [1] at high temperature (T > T c ~ 200MeV), 
which is called as the QCD phase transition. The experimental creation of the 
QGP, one of the most central subjects in the RHIC project, is expected to be re- 
alized in ultra-relativistic heavy-ion collisions [1,25]. On the other hand, the QCD 
phase transition occurred in the early Universe, and its process strongly influenced 
the afterward nucleosynthesis [26]. 

In the DGL theory, the QCD phase transition is characterized by QCD-monopole| 

condensate x = \Xa\, which is an order parameter on the confinement strength. 
The QCD-monopole condensed vacuum (x ^ 0), where the dual gauge symmetry 
is spontaneously broken (mg ^ 0), is the confinement phase with a non- vanishing 
string tension (k ^ 0). Without QCD-monopole condensation (x = 0), the dual 
gauge symmetry is manifest (m# = 0), and the system corresponds to the decon- 
finement phase, where the string tension disappears (k — 0). 

To study the QCD phase transition, the effective potential V e f[(x;T) at finite 
temperature, which physically means the thermodynamical potential, is formulated 
as the function of QCD-monopole condensate x = \Xa\ [10-12,14], 

V eS (x;T) = 3\(x 2 -v 2 ) 2 + 3T J A[2k(l-e-^)da(l-e-^)] 

(4.1) 

with rriB = V3gx an d rn x = 2\f\x- Here, we have used the quadratic source term 
to avoid the imaginary scalar-mass problem [10-12,14]. 

The QCD-monopole condensate Xphys(^) & t finite temperature is obtained from 
the local minimum of the effective potential V e s(x;T). As temperature increases, 
Xphys(T) decreases and disappears at a critical temperature T c . Above T c , the 
deconfinement phase is realized as the Coulomb phase with Xphys 

(T) = 0, where 

the dual gauge symmetry is restored. With the parameters in Ref.[8], we find a 
weak first order phase transition at T c =0.2GeV, and the mixed phase or the two- 
phase coexistence is allowed only in T\ ow < T < T up with T\ ow ~ 0.193GeV and 
T up ~ 0.201GeV [10,12,14]. 

Now, let us consider the confinement and hadron properties. Fig. 4 shows the 
string tension k(T) at finite temperature. As the temperature increases, k(T) 
becomes smaller and drops rapidly near T c [10-12,14,27]. Therefore, the slope of 
the inter-quark potential is reduced and the inter-quark distance inside hadrons 



increases at high temperature. In addition, the color-electric field spreads according 
to the decrease of tub- Thus, the reduction of the confinement force leads to the 
swelling of hadrons at high temperature [14]. 

We predict also a large mass reduction of the scalar glueball originated from 
QCD-monopoles near T c [10-12,14]. We guess that the violent excitation of these 
scalar glueballs with a reduced mass would promote the QCD phase transition. 

Finally, we consider the relation between the surface tension a and the effective 
potential V e g(x',T c ) [10,11]. The surface tension a characterizes the strength of 
the first order in the phase transition, and is very important for the shape of the 
boundary surface in the mixed phase, where the two phases correspond to the two 
minima, x — and x — Xc in Ksff (x T c ). Using the sine-Gordon kink ansatz [10,11] 
for the boundary profile, x(z) = Xctan -1 e z / s , we derive simple formulae for the 

surface tension a ~ -^V^Xc, and the phase-boundary thickness 25 ~ ^^Xc/ V~h, 
where h is the "barrier height" between the two minima in V e ff(x;T c ). We find 
25 ~ 3.4fm and a ~ (112MeV) 3 [10], which seems consistent with the lattice QCD 
data [28]. 

5. Hadron Bubble Formation in the Early Universe 

In this chapter, we study the hadron bubble formation [29] in the early Uni- 
verse using the DGL theory [14]. As Witten pointed out [30], if the QCD phase 
transition is of the first order, the hadron and QGP phases should coexist in the 
early Universe. During the mixed-phase period, there appears the inhomogeneity 
on the baryon density distribution [26,31], which can strongly affects the primordial 
nucleo-synthesis and the successive history of the Universe. 

Let us consider how hadron bubbles appear in the QGP phase near the critical 
temperature T c in the DGL theory. The hadron bubbles are created in the super- 
cooling QGP phase with T[ ow < T < T c . We use the sine-Gordon kink ansatz [14] 
for the profile of the QCD-monopole condensate in the hadron bubble as a function 
of the radial coordinate r, x(r\ R) = Xphys(^) tan -1 e^ R ~ r ^ s / tan -1 e R ^ s , where R 
and 25 correspond to the hadron-bubble radius and the phase-boundary thickness, 
respectively. The QCD-monopole condensate x(r; R) is finite only inside the bub- 
ble, r < R. The energy density £(r;R) of the hadron bubble is shown in Fig. 5. 
It is negative inside and positive near the boundary surface. The total energy is 
roughly estimated as the sum of the positive surface term and the negative volume 
term. 

The total energy of the hadron bubble with radius R can be estimated using 
the effective potential V eS (x;T), E(R; T) = 4tt / °° drr 2 {3(^^) f + V eS (x;T)}, 
where the thickness 5 is determined by the energy minimum condition. The hadron- 
bubble energy E(R; T) as shown in Fig.6 takes a maximal value at a critical radius 
R c , which leads the collapse of the hadron bubbles with smaller radius than R c . 
Only large hadron bubbles with R > R c grow up with radiating the shock wave 
[14,29]. 



On the other hand, the creation of large bubbles is suppressed because of 
the small creation probability [14]. In the hadron-bubble formation, there is a 
penetration over a large barrier height h(T) in the effective potential per unit 
volume, and therefore the creation of large bubbles needs a large energy fluctua- 
tion. Such a process is strongly suppressed because of the thermodynamical factor 
P{T) = exp(— ^-R c (T) 3 h(T)/T), which is proportional to the hadron-bubble for- 
mation rate. Thus, the only small bubbles are created practically, although its 
radius should be larger than R c [14]. 

Using V e s(x;T) in the DGL theory, we can estimate the critical radius R C (T) 
and the hadron-bubble formation factor P(T) at finite temperature as shown in 
Fig. 7 and Fig.8, respectively. As the temperature decreases, the created hadron 
bubbles becomes smaller, while the bubble formation rate becomes larger [14]. 

From these results, we can imagine how the QCD phase transition happens in 
the big bang scenario [14] as shown in Fig.9. 

(a) Slightly below T c , only large hadron bubbles appear, but the creation rate is 
quite small. 

(b) As temperature is lowered by the expansion of the Universe, smaller bub- 
bles are created with much formation rate. During this process, the created 
hadron bubbles expand with radiating shock wave, which reheats the QGP 
phase around them [29]. 

(c) Near T\ ow , many small hadron bubbles are violently created in the unaffected 
region free from the shock wave. 

(d) The QGP phase pressured by the hadron phase is isolated as high-density 
QGP bubbles [29,30], which provide the baryon density fluctuation [26]. 

Thus, the numerical simulation using the DGL theory would tell how the hadron 
bubbles appear and evolve quantitatively in the early Universe. 

We would like to thank Prof. Kajino for fruitful discussions on the hadron 
bubble formation. 
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Table 1: Features of the lattice QCD and the QCD effective models. 

Table 2: The comparison among the ordinary electromagnetic system, the super- 
conductor and the nonperturbative QCD vacuum regarded as the QCD-monopole 
condensed system. 

Fig.l: A brief sketch of the QCD physics. 

Fig. 2: A brief history of QCD and typical QCD effective models. 

Fig. 3 : The electric and magnetic inter-charge potentials in the Coulomb, Higgs 
and dual Higgs phases. The "see-saw relation" is found between in the electric and 
magnetic sectors. The long-range confinement system in one sector corresponds to 
a short-range (Yukawa-type) interaction system in the other sector. 

Fig. 4: The string tension k(T) at finite tem- 
perature T. The black dots denote lattice 
QCD data [27] in the pure gauge. 

Fig. 5: The energy density £(r;R) of the 
hadron bubble. It is negative inside and 
positive near the boundary surface. 

Fig. 6: The energy of the hadron bubble 
E(R] T) as function of the hadron bubble 
radius R. 

Fig. 7: The critical radius R C (T), which pro- 
vides maximum of E(R;T), at finite tem- 
perature T. 

Fig. 8: The hadron-bubble formation factor 
P(T) = exp(-^Rl(T)h(T)/T) as a func- 
tion of temperature T. 

Fig. 9: The appearance and evolution of hadron bubbles in the early Universe. The shaded and 
white regions denote the hadron and QGP phases, respectively, (a) Slightly below T c , only large 
hadron bubbles appear with very small creation rate, (b) Hadron bubbles expand with radiating 
shock wave, which reheats the QGP phase around them. The affected region is expressed by the 
arrow, (c) Near T\ ow , many small hadron bubbles are violently created in the unaffected region, 
(d) The QGP phase pressured by the hadron phase is isolated as high-density QGP bubbles. 



